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Introduction 


Significant  opportunity  to  improve  vehicle 
occupant  safety 

Reduce  impact  severity  between  occupant 
heads  and  vehicle  interiors 

Rigid  body  panels 

Used  plastic  deformation  of  mild  steel  fins 
and  cover  sheet  to  absorb  impact  energy 


Absorber  construction 


•  0.5  inch  wide  mild 
steel  fins 

•  Connected  with  a  mild 
steel  web 

•  Sandwiched  between 
mild  steel  surface 
panel  and  rigid  armor 


LS-DYNA  absorber  model 


Connected  fin  /  web  assembly  to  cover 
sheet  using  spot  welds 

Used  SPC  to  anchor  fins  at  interface 
between  absorber  and  rigid  panel 

Type  13  contacts  between  the  various 
assembly  components 

Nominal  3  mm  mesh 

Steel  was  modeled  using  MAT24 


Impact  test  simulation 


FMVSS  201 U 

Component  level,  10 
inch  X  10  inch  surface 

FTSS  v.  3.6  free 
motion  headform 

15  mph  initial  velocity 

20°  angle  between 
velocity  and  surface 


Optimization  problem 


Minimize  crush  space  subject  to  the 
constraint  that  HIC(d)  <  700 

Independent  variables: 

-  crush  space 

-  spacing  between  fins 

-  shell  thickness  of  fin,  web,  and  cover  sheet 


Head  Injury  Criterion  (HIC) 


HIC  -  max  < 

l\  Ji 


HIC(d)  =  0.75446  (HIC)+ 1 66.4 


•  HIC  is  used  to  estimate  the  severity  of  head  impact  events 

•  HIC(d)  is  a  correlation  between  free  motion  headform  HIC  and  HIC  for 
a  full  50th  percentile  dummy 

•  In  the  expression  for  HIC,  a(t)  is  defined  as  the  resultant  acceleration 
as  a  function  of  time,  t1  and  t2  are  any  two  points  in  time  during  the 
impact  separated  by  not  more  than  36  milliseconds. 

•  Lower  HIC  is  better,  FMVSS  201 U  requires  that  HIC(d)  be  less  than 
1000 


Optimization  technique 

Closely  followed  Stander  and  Craig  (LS-OPT) 
successive  response  surface  method 

Iterative  sequence  of  linear  least  squares 
response  surfaces 

Chose  D-Optimal  subsets  of  33  full  factorial 
basis  designs 

D-Optimal  subsets  contained  seven 
combinations  of  the  three  design  factors 

15  iterations  of  7  runs  each;  105  simulations 
overall 


Convergence  to  optimum  values 


Initial 

Optimum 

Crush  space  (inch) 

0.875 

0.8044 

Fin  spacing  (inch) 

0.875 

0.9446 

Fin  /  web  /  cover  shell  thickness  (inch) 

0.0285 

0.02616 

HIC(d) 

737 

699 

Convergence  -  crush  space  and 

spacing 


Convergence  -  shell  thickness 


Acceleration  -  time  history  for 

optimum  design 


Variability  /  robustness 


Optimum  design  settings  ±  5% 


Initial  design  interaction 


Crush  Space 


Interaction  Plot  (data  means)  for  HIC(d) 


0.8974  0.9446  0.9918  0.02485  0.02616  0.02747 


Crush 

Space  Point  Type 
0.7642  Corner 

0.8044  Center 

0.8446  Corner 

Spacing  Point  Type 
0.8974  Corner 

0.9446  Center 

0.9918  Corner 

Thickness 


Crush  space  -  thickness  interaction 


t  (s) 


low  crush  space  -  low  th  HIC(d)=818 
low  crush  space  -  hi  th  HIC(d)=756 
hi  crush  space  -  low  th  HIC(d)=709 
hi  crush  space  -  hi  th  HIC(d)=688 


•  Lower  shell  thickness  increases  propensity  toward 
“bottoming  out” 

•  Lower  crush  space  tends  to  increase  mean  deceleration 


HIC(d)  =  688,  first  peak 


SAE  J1000  Filtered  Acceleration 


HIC(d)  =  688,  second  peak 


0  0  002  0  00*  0  0C5  0  008  0  01  0  012  0  014  0.01©  0  018  0  02 

Time  (s) 


Improve  design 


Cube  Plot  (data  means)  for  HIC(d) 


©  Centerpoint 
•  Factorial  Point 


-5%  increase  in  crush  space  should  enable  HIC(d)  <  700 
Use  0.84  crush  space  face  as  starting  point 


Design  improvement 


Used  interpolation  process  keeping  spacing  and  thickness  fixed 

Found  new  crush  space  to  ensure  HIC(d)  <  700  for  nominal 
parameter  settings  ±  5% 


Design  improvement 


Improve  design 


Cube  Plot  (data  means)  for  HIC(d) 


©  Centerpoint 
•  Factorial  Point 


-5%  increase  in  crush  space  should  enable  HIC(d)  <  700 
Use  0.84  crush  space  face  as  starting  point 


Improved  design 


Cube  Plot  (data  means)  for  HIC(d) 


©  Centerpoint 
•  Factorial  Point 


New  nominal  design  settings  ±  5% 

Moderate  (~  0.1  inch)  increase  in  nominal  crush  space  yields 
HIC(d)  <  700 


Improved  design  interaction 


Crush  Space 


Interaction  Plot  (data  means)  for  HIC(d) 


0.8974  0.9446  0.9918  0.02485  0.02616  0.02747 


Crush  Space  Point  Type 
0.85262  Corner 
0.89749  Center 
0.94236  Corner 


Spacing  Point  Type 
0.8974  Corner 
0.9446  Center 
0.9918  Corner 


Thickness 


Response  surfaces  for  the 
improved  design 

Sampled  by  means  of  uniform  designs 

Developed  response  surfaces  via  Kriging 

Factorial  simulation  results  were  used  to 
compare  fidelity  of  Kriging  response 
surfaces  generated  in  various  ways 


Kriging 


f„(x)  +  Z(x) 


R  (  vf  -  (,  I  =  exp  [  -  ^  I  a,  -  a  . )’  ] 


Kriging  models 

Compared  results  for  surfaces  generated 
with 

-  constant 

-  first  order  polynomial 

-  quadratic  polynomial 

Gaussian  correlation  function 

Three  different  sample  sizes  -  9,  17,  and 
30 


Goodness-of-fit  estimates 


Maximum  Error  =  max 


yk 


nge ,  i 


y factorial. 


RMSE  = 


y  factorial 


1/2 


Comparison  of  maximum  error 


Sample  size 

Constant 

First  Order 
Polynomial 

Quadratic 

Polynomial 

9 

42.70 

29.52 

17 

125.28 

37.44 

34.92 

30 

53.39 

30.08 

18.66 

Comparison  of  root  mean  square 

error  (RMSE) 


Sample  size 

Constant 

First  Order 
Polynomial 

Quadratic 

Polynomial 

9 

24.36 

18.36 

17 

60.50 

16.00 

19.02 

30 

23.95 

15.04 

10.67 

Kriging  model  contours 


Spacing  =  0.95  (5%  lower  than  nominal  value) 
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Kriging  model  contours 


Spacing  =  1.00  (nominal  value) 
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Kriging  model  contours 


Spacing  =  1.05  (5%  higher  than  nominal  value) 
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Conclusions 

It’s  possible  to  very  efficiently  optimize  an  energy 
absorber  design  using 

-  the  LS-DYNA  explicit  finite  element  code 

-  the  successive  response  surface  method  algorithm 

Use  of  classic  factorial  techniques  in  combination  with 
Kriging  response  surfaces  can 

-  guide  improvement  of  product  robustness 

-  offer  insight  into  the  nature  of  a  product  and  its  performance 
variability 

An  enlightened  combination  of  these  techniques 
enables,  if  nothing  else,  valuable  and  relatively 
inexpensive  insight  into  the  feasibility  and  behavior  of 
various  design  concepts. 


